Feedback spin resonance in superconducting CeCu2Si2 and CeCoIns 
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We show that the recently observed spin resonance modes in heavy- fermion superconductors 
CeCoIns and CeCu2Si2 are magnetic excitons originating from superconducting quasiparticles. The 
wave vector Q of the resonance state leads to a powerful criterion for the symmetry and node posi- 
tions of the unconventional gap function. The detailed analysis of the superconducting feedback on 
magnetic excitations reveals that the symmetry of the superconducting gap corresponds to a singlet 
d^2_y2 state symmetry in both compounds. In particular this resolves the long-standing ambiguity 
of the gap symmetry in CeCoIns . We demonstrate that in both superconductors the resonance peak 
shows a significant dispersion away from Q that should be observed experimentally. Our results 
suggest a unifying nature of the resonance peaks in the two heavy-fermion superconductors and in 
layered cuprates. 

PACS numbers: 74.20.-z, 74.70.-b 



The relation between unconventional superconductiv- 
ity and magnetism in heavy-fermion systems and doped 
transition metal oxides is one of the most interesting re- 
search areas in condensed matter physics. Despite of 
certain differences concerning the proximity to a Mott 
insulator in the transition metal oxides and the weak hy- 
bridization of the /-electrons in the heavy-fermion sys- 
tems, it is widely believed that in both cases short-range 
antiferromagnetic (AF) spin fluctuation are responsible 
for Cooper-pairing with a d-wave order parameter. Fur- 
thermore, unconventional superconductivity itself has a 
strong feedback on the magnetic spin excitations in these 
systems below the superconducting transition tempera- 
ture Tc. One example is the famous resonance peak ob- 
served in high-Tc cupmtes by means of inelastic neu- 
tron scattering(INS) [l] whose nature is still actively 
debated [2]. The complexities of the cuprates including 
the electronic inhomogeneities or stripe formation com- 
plicates here the interpretation of the experimental data. 

For some time the cuprates were considered to be 
unique in displaying the resonance peak in the supercon- 
ducting state. Then, remarkably INS revealed the forma- 
tion of a new magnetic mode in the superconducting state 
of the 5/-heavy fermion compound UPd2Al3 with Tc = 
1.8 K y. Its sharply peaked intensity, its temperature 
dependence and energy position well below 2Ao (with Aq 
being the maximum of the superconducting gap) strongly 
resembles the resonance peak seen in high-Tc cuprates. 
In distinction to the cuprates, the symmetry of the gap 
function was unclear for UPd2Al3. Thermal conductivity 
measurements in a rotating magnetic field showed that 
it has nodes in the hexagonal plane but could not de- 



As already noticed 
the observation of a 



termine the symmetry uniquely 
in [j, l5[ and shown explicitly in 
resonance at the AF wavevector Q puts a stringent con- 
dition on the gap symmetry by requiring a sign change 



Ak+Q = — Ak under translation by Q. Together with 
thermal conductivity results this determines unambigu- 
ously the gap function as Ak = Aq cos ckz which has node 
lines at the boundaries of the antiferromagnetic Brillouin 
zone. The case of UPd2 AI3 shows not only that the mag- 
netic resonance due to a feedback effect of the supercon- 
ducting state is a universal phenomenon in unconven- 
tional superconductors but also, that inelastic neutron 
scattering (INS) is a powerful technique to determine the 
superconducting gap symmetry. 

Indeed, very recently the resonance peak has been also 
observed below the superconducting transition temper- 
atures in Ce-based heavy-fermion compounds, namely 
in CeCu2Si2[7] at the (incommensurate) wave vector 
QsDW ^ (0.22 • 27r/a,0.22 • 27r/a,0.53 • 27r/c) and in 
CeCoIns [8] at the antiferromagnetic wave vector Qaf ^ 
(yr/a, 7r/a, tt/c) . The latter system which has quasi-two- 
dimensional tetragonal crystal structure shows the high- 
est superconducting transition temperature Tc = 2. IK 
among heavy-fermion compounds [9|. Its gap symmetry 
has been long disputed because of conflicting results from 
angle resolved magnetothermal conductivity [lO| {dj.2_y2 
-symmetry) and specific heat [Uj (c^^^^-symmetry) mea- 
surements. The INS results [8[ and the analysis presented 
here give a clear resolution of this puzzle in favor of the 
dx2_y2 state which underlines again the importance of 
INS for determining the superconducting gap symme- 
try. Heavy quasiparticles in the Ce compounds [l2|, ll3[ 
of predominantly 4/ character are more strongly corre- 
lated than 5/ quasiparticles or quasiparticles in cuprates. 
Therefore, similarities of the resonance state indicate a 
generic feature. 

Here, we analyze the dynamical magnetic susceptibil- 
ity in CeCu2Si2 and CeCoIns below the superconducting 
transition temperature, assuming unconventional charac- 
ter of the Cooper-pairing originating from the exchange 



of the antiferromagnetic spin fluctuations. We show that 
in both cases the resonance feature evolves at the mag- 
netic instability wave vector. We discuss the dispersion of 
the resonant excitations as a function of the momentum 
and demonstrate that in both compounds the dx2_y2- 
wave symmetry of the superconducting order parameter 
is consistent with the INS experiments. 

Starting our analysis by considering the CeCu2Si2 sys- 
tem, we note that a good starting point for calculat- 
ing the magnetic susceptibility is the band structure ob- 
tained within the renormalized band theory [ij, [isl, llq . 
In this approach the phase shifts for the 4/-states are in- 
troduced empirically into the ab-initio local density ap- 
proximation (LDA) scheme to generate heavy quasipar- 
ticle bands at the Fermi level. Crystalline Electric Field 
(CEF) splitting strongly affects the quasiparticle disper- 
sion and has to be properly accounted for [l7[. Note 
that the formation of a resonance in the spin response 
is determined mainly by the unconventional symmetry 
of the superconducting order parameter. Therefore it is 
reasonable to use a tight-binding parametrization fit to 
the main heavy quasiparticle band that crosses the Fermi 
level. The corresponding Fermi surface consisting of the 
stacked columns along c direction is shown in FigUJa) in 
the first few Brillouin zones. The obtained Fermi surface 
shows a fiat part connected by the nesting wave vector 
QsDW = (0.22 • 27r/a, 0.22 • 27r/a, 0.52 • 27r/c) as indicated 
by the arrow. As discussed in earlier work [18| the Qsdw 
wavevector agrees very well with the experimentally ob- 
served SDW in the so-called A-phase of that compound. 



(b) q.=0.52 r.l.u 
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FIG. 1: (a) Calculated Fermi surface for the main electron 
sheet in CeCu2Si2 using the tight-binding parametrization in 
the form £k = 2ti (cos kxa + cos kya) + 4t2 cos kxa cos kya + 
8t3 cos ^ cos %^ cos -^ -/i where ti = 17.5, ts = -5.2, 
ts = —11.2 and /a = —57.4 (in K) are the hopping integrals 
and the chemical potential, respectively. Here, a = 4.1 A and 
c = 9. 92 A are the lattice constants. The arrow indicates 
the SDW scattering wave vector, (b) Calculated static spin 
susceptibility on two-dimensional mesh for qz = 0.52 r.l.u. 
The arrow indicates the SDW ordering wave vector, Qsdw, 
as observed in the experiment J14l ]. 

In Figlljb) we show the real part of the Lindhard 
susceptibility Xo(q) calculated with the parameterized 
heavy quasiparticles. One finds that Re x{^^ 0) is peaked 
at Qsdw- When compared with a fully renormalized 



band structure calculation [1^ we observe that the tight- 
binding bands result in a somewhat more pronounced 
nesting and yield a whole contour of the nesting wave 
vectors around Qsdw- This is due to the relatively sim- 
ple band structure that contains nearest neighbor hop- 
ping integrals only. However, this does not infiuence our 
results concerning the resonance feature below Tc, since 
the latter is determined mainly by the special nature of 
the superconducting gap. 

The resonance peak in the SC state of CeCu2Si2 as 
well as of CeCoIns can be understood by considering the 
dynamical spin susceptibility within the random phase 
approximation (RPA), ie.. 



Xi?pA(q,^) = 



xo(q,^) 



i-^qXo(q,^)' 



(1) 



where 17^ is the fermionic four-point vertex and Xo(q, ^) 
is the heavy quasiparticle susceptibility. The latter is 
given by the sum of the well-known bubble diagram con- 
sisting of either normal or anomalous (T<Tc) Green 
functions. For large momenta q, Imxo(q, 0) is zero at 
low frequencies and can exhibit a discontinuous jump 
at the onset frequency of the p-h continuum Qc = 
min(|Ak| + lAk+al) where both k and k+q lie on the 
Fermi surface[l9|. Note, however, that the discontinuity 
in Imxo occurs only if sgn(Ak) = — sgn(Ak+q) which is 
only possible for unconventional order parameters. A 
discontinuity in Imxo leads to a logarithmic singular- 
ity in Rexo- As a result, the resonance conditions (i) 
/7qRexo(q,^res) = 1 and (ii) Imxo(q,^res) = can be 
both fulfilled at uOres < ^c for any /7q > 0, leading to the 
occurrence of a resonance peak in form of a spin exciton 
below Tc. For finite quasiparticle damping F, condition 
(i) can only be satisfied if /7q > exceeds a critical value, 
while condition (ii) is replaced by ImxoCq, ^res) << 1. 

Since the symmetry of the superconducting gap has 
not yet been determined unambiguously in CeCu2Si2, we 
have analyzed all the spin singlet s- and d-wave func- 
tions allowed by the crystal-group symmetry of the lat- 
tice [IJ, [20, l2l|. We have found that the resonance- 
like feature, i.e., the discontinuous jump in Imxo and 
the corresponding logarithmic singularity in Rexo oc- 
cur at the SDW wave vector Qsdw for three types 
of the order parameters: Ak = Aq (cos kxa — cos kya) 
belonging to the Big irreducible representation, Ak = 
Ao sin ^ sin ^ cos ^ belonging to the B2g irreducible 
representation, and also for each of the two components 
of the Eig representation, Ak = Aq sin/ca^asin/c^c and 
Ak = Aq sin ^ "" ^y^^ sin ^. An important finding is 
that for a given dispersion the resonance in the Big 
channel is by far the strongest. In Fig|2](a) we show 
the RPA susceptibility in the normal and the supercon- 
ducting state for the Big (Ak = Aq {cos kxa — cos kya)) 
and for the B2g (Ak = Aq sin ^ sin -^) channels, re- 
spectively. One finds that the susceptibility for the Big 



symmetry is larger in the superconducting state than it 
is in the normal state, while in the B2g channel there 
is no enhancement of the normal state spin susceptibil- 
ity. This definitely points towards a dx2_y2-wdive sym- 
metry of the superconducting gap in CeCu2Si2 since a 
very sharp resonance was found at QsDW in INS [7|. We 
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FIG. 2: (a) Calculated real and imaginary parts of the RPA 
spin susceptibility for the normal (red) and superconducting 
state for the B\g (blue) (Ak = Ao (cos /ccca — cos /c^a)) and 
B^g (green) (Ak = Ao sin ^|^ sin ^^) symmetry of the su- 
perconducting gap. Here, we assume Uqisdw ~ ^ti to satisfy 
the resonance condition in the superconducting state. For 
the numerical purpose we also set the damping F = 2K. (b) 
Calculated dispersion of the resonance peak in CeCu2Si2 for 
(0.22, 0.22, qz) direction. The dispersion is nearly flat around 
^SDW due to the two-dimensional structure of the supercon- 
ducting gap. Here, we have introduced an interaction peaked 



at QsDW in the form [/q 
with B=0.5. 
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note, however, that the magnitude of the resonance for 
each symmetry is a result of a competition of various 
features, i.e., the curvature of the Fermi surface at the 
points connected by Qs-dvf, the absolute value of the 
superconducting gap and also the velocity at the node 
of the gap. Therefore a modification of the electronic 
dispersion, and the actual angular dependence of the su- 
perconducting gap can modify the results for the absolute 
intensity of the resonance peak. This is, for example, the 
case for the non-monotonic d— wave gap in the electron- 
doped cuprates[22|. 

The resonance condition is also satisfied for momenta 
Q^ slightly away from ^sdw as long as sgn(Ak) = 
— sgn(Ak+Q.). We find that the resonance is readily sup- 
pressed by variation of the in-plane (g^^, g^y) -momentum 
and exists only in the close vicinity to ^SDW- This 
is because for an incommensurate momentum there are 
always scattering processes which involve parts of the 
Fermi surface with sgn(Ak) = +sgn(Ak+Q.) thus sup- 
pressing the resonance. In addition, varying Q^ away 
from (^SDW along the in-plane (qx^qy^Qz^^) direction 
the nesting condition is also lost which overall yields a 
decrease of Rexo (see also Figdjb)). Since the suppres- 
sion of the resonance occurs for all symmetries of the 
above mentioned superconducting order parameters the 
resonance peak is confined to momentum {qx^Qy^Qz^^) 
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FIG. 3: (a) Calculated Fermi surface for CeCoIns using the 
band structure parameters adopted previously [23^1 . The dash- 
dotted arrow points at states at the Fermi surface scattered 
by the antiferromagnetic wave vector, Qaf. The dashed lines 
depict the position of nodes in the first BZ for a superconduct- 
ing order parameter of d^2_y2-waYe symmetry . Following 
Ref. [23 we set the energy unit 0.26eV. (b) Calculated static 
spin susceptibility of non-interacting electrons in CeCoIns for 
qz = 0.5 r.l.u. on a two-dimensional mesh. 



in the plane around Qsdw- At the same time the situ- 
ation is less dramatic with respect to the qz momentum 
dependence. In Fig|2](b) we show the dispersion of the 
resonance along the qz direction. One finds that it re- 
mains nearly fiat as one departs from Qsdw- As a mat- 
ter of fact, for a constant interaction /7q, the resonance 
shows a weak dispersion resulting from the slight change 
of the Rexo away from Qsdw- 

Let us now turn to the resonance peak formation 
in CeCoIn5 [pi . According to the band structure 
calculations [23|, CeCoIns comprises several /- and con- 
duction bands which are hybridized in a complex manner. 
Therefore, it is quite difficult to reproduce a resulting 
structure by using a single band model. Introducing a 
two-band model is much more appropriate. In fact, this 
has been done previously for CeCoIns [23] and it has 
been found that the resulting energy dispersion crossing 
the Fermi level can be written as: 



^2k = - 



(4 + ^^)-V(^k^-£ 



c)^- 



w^ 



(2) 



where E^ and e^ is the effective /—band and the conduc- 
tion band dispersions, respectively, and Vjk is the effec- 
tive hybridization strength, renormalized by the on-site 
/ — /-Coulomb repulsion. The resulting Fermi surface is 
shown in Fig. [3ja). Like for CeCu2Si2 the present Fermi 
surface has again nesting properties. However, here it 
occurs for the the commensurate antiferromagnetic wave 
vector Qaf = (7r/a,7r/a,7r/c). This agrees with recent 
INS data on the normal state [8(. In FiglSjb) we show 
the real part of the Lindhard spin susceptibility for the 
normal state calculated on a two dimensional mesh. In 
accordance with the Fermi surface topology we find that 
the spin susceptibility is peaked at the antiferromagnetic 
wave vector, Qaf- R is important that these antiferro- 
magnetic spin fluctuations are responsible for the forma- 
tion of d-wave superconductivity in this compound. 



To address the issue of the resonance peak forma- 
tion, we show in FiglJl^a) the calculated real and imagi- 
nary part of the RPA susceptibility at the antiferromag- 
netic wave vector in the normal and the superconduct- 
ing states. Among possible superconducting symmetries 
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FIG. 4: (a) Calculated real and imaginary part of the RPA 
susceptibility at the antiferromagnetic wave vector Qaf for 
CeCoIns as a function of frequency in the normal (red) and 
superconducting, Bi^ (blue) and B2g (green) states. We have 
assumed that t/QAF = f^o ~ 1.66t. (b) Calculated disper- 
sion of the resonance peak in CeCoIns for the d^2_y2-waYe 
symmetry of the superconducting order parameter along the 
{Qx^qy^Tv/c) direction. The abscissa is shown in units of vr/a. 
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in CeCoIns a resonance peak forms only for the Big ( 
Ak = -Y-{coskxa — cos ky a)) symmetry of the supercon- 
ducting order parameter. As in the case of CeCu2Si2 
the antiferromagnetic wave vector connects states with 
opposite sign of the superconducting gap. As is clearly 
visible from Figl3l(a) this results in the formation of a 
resonance peak similar to the one in CeCu2Si2. The 
resonance peak forms near the particle-hole continuum, 
i.e., close to Qc = min(|Ak| + |Ak+q|) which is around 
Aq. This is because the points connected by Qaf are 
lying relatively far from the part of the Fermi surface 
where the gap function has maximum value. Note, if 
the Cooper-pairing itself arises due to an exchange of an- 
tiferromagnetic spin fluctuations, the maximum of the 
superconducting gap occurs at points of the Fermi sur- 
face which are connected by Qaf- At the same time, the 
symmetry of the superconducting gap possesses still the 
(i^2_^2-wave symmetry, although with higher harmonics 
included. It is remarkable that like in CeCu2Si2 we find 
that only a gap function of dx2_y2-wdive (Big) type re- 
sults in the formation of the resonance peak at Qaf- 
This unambiguously confirms the bulk symmetry of the 
superconducting gap in CeCoIns. The dxy-wdive {B2g) 
symmetry discussed in the literature is clearly ruled out. 
Finally, in Fig. [41(b) we show the dispersion of the 
resonance excitations away from the Qaf- We observe 
that the resonance disperses downwards as a function 
of frequency. This behavior is similar to that found in 
hole-doped high-Tc cuprates. In particular, the value 
of the critical frequency Qc = |Ak| + |Ak+q| lowers for 
q < Qaf for the in-plane momentum, since the scatter- 



ing occurs for states closer to the diagonal part of the 
Brillouin zone where the superconducting gap is smaller. 
As a result, the resonance condition shifts to lower ener- 
gies. Our results show that the resonance peak possesses 
a universal dispersion away from the wave vector that 
connects the maxima of the superconducting gaps and 
we suggest to investigate experimentally the dispersion 
of the resonance. 

In conclusion, we analyze the dynamical magnetic sus- 
ceptibility in CeCu2Si2 and CeCoIns below supercon- 
ducting transition temperature. We show that in both 
cases a resonance feature evolves at the wave vector of 
the magnetic instability. Our results show that the lat- 
ter two heavy-fermion superconductors and the high-Tc 
cuprates possess the same symmetry of the supercon- 
ducting order parameter suggesting that the same mecha- 
nism of the Cooper-pairing is probably involved. Further- 
more, despite of the three-dimensional electronic struc- 
ture, the two-dimensional (ia,2_^2-wave superconducting 
gap in CeCu2Si2 and CeCoIns may provide further hints 
on the microscopic mechanism of unconventional super- 
conductivity in heavy-fermion systems as well as in lay- 
ered cuprates. 
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